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We have used the Space Telescope Imaging Spectrograph (STIS) aboard HST to measure 
a sulfur abundance of [S/H] = -0.41 ± 0.06 in the interstellar medium (ISM) of the nearby 
damped Lyman-a (DLA) absorbing galaxy SBS 1543+593. A direct comparison between this 
O , QSO absorption line abundance on the one hand, and abundances measured from H II region 

emission line diagnostics on the other, yield the same result: the abundance of sulfur in the neu- 
| tral ISM is in good agreement with that of oxygen measured in an H II region 3 kpc away. Our 

result contrasts with those of other recent studies which have claimed order-of-magnitude differ- 
ences between H I (absorption) and H II (emission) region abundances. We also derive a nickel 
abundance of [Ni/H] < -0.81, some three times less than that of sulfur, and suggest that the de- 
pletion is due to dust, although we cannot rule out an over-abundance of alpha-elements as the 
cause of the lower metallicity. It is possible that our measure of [S/H] is over-estimated if some 
S II arises in ionized gas; adopting a plausible star formation rate for the galaxy along the line of 
sight, and a measurement of the C II* A 1335. 7 absorption line detected from SBS 1543+593, we 
determine that the metallicity is unlikely to be smaller than we derive by more than 0.25 dex. We 
estimate that the cooling rate of the cool neutral medium is log [l c (ergs s -1 H atom -1 )] rs -27.0, 
the same value as that seen in the high redshift DLA population. 

Subject headings: quasars:absorption lines — quasars individual (HS 1543+5921) 
— galaxies:individual(SBS 1543+593) — galaxies:abundances 



1. Introduction 



Since the launch of the Hubble Space Telescope (HST) it has been hoped that a comparison could be 
made between galactic elemental abundances deduced, respectively, from emission line measurements made 



'Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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of H II regions within some designated nearby galaxy, and ultraviolet (UV) absorption line data recorded 
in the spectra of a UV-bright source whose sight line passed through the same galaxy. Both measurements 
would yield the composition of the present-day interstellar medium (ISM) of such a galaxy, either in H II 
regions ionized by OB stars, or in neutral, diffuse, interstellar clouds which happen to lie in front of the 
background source. The initial expectation has been that the two techniques should give the same answer 
when applied to the same location within a particular galaxy. Indeed, for the solar neighborhood the agree- 
ment has now been verified empirically to a satisfactory degree of precision (see the discussion of this point 
by Esteban et al. 2004). Why, then, should we be concerned that this may not be the case in other galaxies? 

One reason arises from the puzzling result that the metallicity of damped Lya systems (DLAs) — the 
class of QSO absorbers with the highest column densities of neutral gas — is consistently below solar at 
essentially all redshifts, from z = 0.1 to z > 4 (Kulkarni et al. 2005; Prochaska et al. 2003; Pettini et al. 
1999). If DLAs were representative of the galaxy population as a whole (e.g., Wolfe 1990), their element 
abundances should grow with time to the near-solar levels typical of today's galaxies (Fukugita & Peebles 
2004). Radial abundance gradients have been appealed to as a possible explanation for the lower than solar 
values (Christensen et al. 2005; Chen et al. 2005; Ellison et al. 2005). Quantitatively, however, the magnitude 
of such gradients in nearby spirals has recently been questioned (Bresolin et al. 2004), and remains unknown 
at high redshifts. 

A different aspect of the same 'problem' is the possibility that the H II regions of star-forming dwarf 
galaxies may be self-enriched in the nucleosynthetic products — chiefly oxygen — of the most massive, and 
short-lived, stars formed during a bout of star formation. This proposal, originally put forward by Kunth 
& Sargent (1986), has recently been reconsidered in the light of Far Ultraviolet Spectroscopic Explorer 
(FUSE) observations of a number of low metallicity dwarf galaxies, summarized by Cannon et al. (2005). 
The FUSE data seem to show that the abundances of both oxygen and nitrogen in the absorbing H I gas are 
systematically lower, by up to one order of magnitude, than those deduced from nebular emission lines in H II 
regions. Since nitrogen, a product of intermediate mass stars, is released into the ISM several hundred Myr 
after oxygen, such an offset may indicate the existence of an extended, lower-metallicity halo surrounding 
the inner star-forming regions of dwarf galaxies, rather than being evidence of self-pollution by individual 
H II regions. 

The seemingly straightforward consistency check of measuring abundances in the spectrum of a QSO 
lying behind a galaxy with good emission-line metallicity diagnostics has proved deceptively difficult to 
carry out, since the number of suitable QSO-galaxy alignments is small: the foreground galaxy should - 
ideally - be at low redshift so that abundances in individual H II regions can be measured from the ground; 
and the background source must be bright enough in the UV to be observed at high resolution with HST, 
since the important diagnostic lines lie in that wavelength range. Although good QSO-galaxy alignments 
exist, most of the UV sources lie 1 - 2 magnitudes below the limit reachable with HST to make this technique 
routine. 

An alternative method has been to use previously identified intermediate redshift QSO absorption sys- 
tems and to measure the metallicity of the putative absorbing galaxy. In this case, light is collected across 



-3- 



the entire galaxy, rather than from individual H II regions, which are unresolved at such redshifts. Chen 
et al. (2005) conducted such a test in a few galaxies believed to be the hosts of low-redshift DLAs. Again, 
although there are suggestions of an offset, its origin is open to different interpretations, partly because the 
element measured in absorption is iron, which exhibits a highly variable degree of dust depletion in the 
interstellar medium (as well as having a different nucleosynthetic origin from oxygen). 

The metallicity of a galaxy is, of course, closely related to the formation and evolution of its stellar 
population and the consequent enrichment of the interstellar medium within the galaxy. A recently devel- 
oped line of inquiry into the star formation process and the regulation of the interstellar gas in DLAs has 
been through measurements of C II* absorption lines. The ground state of C II is split into two fine-structure 
levels, and the transition from one of these [the (2p) 2 P^/2 level] to the first excited ( 2 D) state produces a 
C II* absorption line at 1335.7 A. The population of the fine structure levels is determined by collisions 
of the C II ion with electrons and neutral hydrogen atoms, and the line can be used to measure important 
diagnostics. Wolfe et al. (2003b) used C II* observed in 30 high-z DLAs to infer a star formation rate (SFR) 
of 6 x 10~ 3 M & yr" 1 kpc" 2 for C II* arising predominantly in a cold neutral medium (Wolfe et al. 2003a, 
2004). They were also able to measure a cooling rate of ~ 10~ 27 ergs s -1 H atom -1 , a value one dex smaller 
than that found by Lehner et al. (2004) for gas cooling in the Milky Way. 

In this paper we present a direct comparison between emission- and absorption-based abundance de- 
terminations in a nearby galaxy and report that, unlike the previous studies reviewed above, we do find 
excellent agreement between the two sets of metallicity measurements. 

In a previous paper (Bowen et al. 2001b) we found that the dwarf galaxy SBS 1543+593 (at a helio- 
centric systemic velocity of 2868 km s -1 — see Bowen et al. 2001a) produces a DLA in the spectrum of 
the background QSO HS 1543+5921 (z em = 0.807; Reimers & Hagen 1998). The QSO sight line intersects 
the inner regions of the galaxy, at only 2.4arcsec, or ~ 0.5 fc^ kpc, from its visual center 1 . Using low- 
resolution spectra from HST, we found that the galaxy produced a DLA with Af(H 1) = 10 2035 cm" 2 . From 
their analysis of the emission line spectrum of the brightest H II region located in a faint spiral arm 3.3 kpc 
from the center of the galaxy, Schulte-Ladbeck et al. (2004) concluded that the abundances of oxygen and 
nitrogen are approximately 1/3 and 1/10 solar respectively, consistent with the low luminosity {Mb = -16.8) 
of SBS 1543+593. In a subsequent paper, Schulte-Ladbeck et al. (2005) determined the sulfur abundance 
of the same emission line region to be ~ 0.5 ± 0.4 times the solar value. They also used the absorption line 
data presented in this paper to compare the metallicity of the neutral gas with that from the H II region, and 
determined that emission and absorption line measurements yield similar abundances. 

In this paper, we present a detailed analysis of our data, taken with the Space Telescope Imaging 
Spectrograph (STIS) aboard HST. We refine our initial estimate of Af(H I), determined originally from lower 
resolution STIS data, and calculate the elemental abundance of sulfur and a lower limit for the abundance of 
nickel. We also investigate possible ranges for the abundances of oxygen, nitrogen and silicon, and highlight 



l h = Ho/70, where Ho is the Hubble constant, and qo = -0.55 is assumed throughout this paper. No corrections are made for 
the velocity of the Milky Way relative to the Cosmic Microwave Background. 
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the uncertainty in the results which arise from using the available STIS data. Of particular concern is the 
question of how much of the observed S II absorption might actually arise in ionized gas, and not just from 
the neutral component measured from the Lya line. We perform conventional photo-ionization modeling 
to show that most of the S II comes from the same neutral component traced by the H I; however, we 
also develop a new method to calculate an ionization correction based on a plausible measure of the star 
formation rate along the line of sight, and the partial measurement of the C II* A1335.7 line found to arise 
in SBS 1543+593. This method again predicts that most of the S II arises in the H I absorbing medium. We 
conclude by comparing our abundances with those measured in the H II regions of SBS 1543+593, as well 
as those found in DLAs. We comment briefly on the fact that the cooling rate of the cool neutral medium in 
SBS 1543+593 may well be similar to that of the DLAs. 

2. Data Acquisition and Analysis 
2.1. Observations 

An ideal observation of HS 1543+5921 would use an echelle mode of STIS to record many interstellar 
absorption lines simultaneously. Unfortunately, the QSO is too faint to be observed in a reasonable allocation 
of HST time. However, HS 1543+5921 is located in the continuous viewing zone (CVZ) for HST, which 
made it possible to observe using the G140M medium-resolution grating. We used the G140M at three 
wavelength settings (see Table 1) to record 55 A- wide portions of the ultraviolet spectrum using the 52x0.1 
aperture. The resolution of these observations is ~ 20 km s" 1 (Kim Quijano 2004). To reduce the effects 
of fixed-pattern noise, data were taken using the dither pattern STIS-ALONG-SLIT, which moved the 
QSO along the slit in increments of 0.55 arcsec. As a result, spectra were recorded at one of three possible 
positions on the detector. 

In choosing which atomic transitions to cover, we were guided by the low resolution (G140L) spectrum 
of the QSO (Bowen et al. 2001b) which showed that the O I A1302 line and the N I AA1199, 1200.2, 1200.7 



Table 1 . Journal of Observations 1 





Exposure Time 


CVZ 


Central Wavelength 




Date 


(s) 


Orbits 


(A) 


Absorption Lines Covered 


2003 Aug 12 


25 250 


5 


1222 


Lya; N I AA1 199, 1200.2, 1200.7; Si III A1206 


2003 Aug 10 


39550 


8 


1272 


S II AA1259, 1253, 1250; Si II A1260 


2004 Jan 29 


14400 


3 


1272 


S II AA1259, 1253, 1250; Si II A1260 


2004 Jan 27 


49500 


10 


1321 


O I A1302; Ni II A1317; C II A1334; 



CH* A1335; Si II A 1304 



All the observations were conducted under GO program 9784. All data were taken with the 52x0. 1 aperture. 
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triplet are likely to be saturated (as is usually the case in DLAs). Furthermore, at zdla = 0.0096, the N I 
lines fall on the wing of the Milky Way damped Lya absorption line, where the flux is reduced to only 
~ 1 x 10" 15 ergs cm" 2 s -1 A" 1 , too low to enable us to record these features with adequate S/N. 

For these reasons, our first tilt setting of the G140M grating was chosen to record the S II AA1259, 
1253, 1250 triplet which we expected to be weak and unsaturated, yet detectable in a reasonable exposure 
time. Sulfur is a good proxy for oxygen; like oxygen (and nitrogen) it is largely undepleted by dust in the 
Galactic ISM and, being an a-capture element, its abundance tracks that of oxygen in Galactic metal-poor 
stars (Nissen et al. 2004). Ionization corrections are expected to be small for both O I and S II when N(H I) 
is large. 

Our second tilt was set to record a relatively short exposure with the grating centered at 1222 A to cover 
the Lya absorption line from SBS 1543+593, and thereby measure a precise H I column density. The N I 
triplet absorption was observed at the same time, and while the low S/N was not expected to be a major 
concern for the damped Lya line (which covers many resolution elements), we knew that the short exposure 
time would make our observations of the N I triplet of little use. 

Finally, we included observations centered at 1321 A to cover the Ni II A1317 absorption line. This 
line is of interest because nickel is a highly depleted element in local interstellar gas. Thus, the Ni II/S II 
ratio, when compared to its solar value, should give an indication of the degree to which refractory elements 
are depleted onto dust in the ISM of SBS 1543+593. This grating tilt also included O I A1302, C II A1334, 
and C II* A1335 absorption from the dwarf galaxy. These lines were expected to be of limited use, but as 
we show below, we were actually able to use them to impose some constraints on physical conditions in the 
absorbing gas. 

2.2. Data extraction 

The pipeline reduction of the spectra centered at 1222 A proved inadequate for the smooth removal of 
the geocoronal Lya emission line, and provided an incorrect estimate of the background — significant counts 
were observed at the bottom of the Lya profiles, where the flux was expected to be zero. Since these prob- 
lems arise from poor extraction of one-dimensional (ID) spectra from the two-dimensional (2D) frames, we 
embarked on an optimal extraction using the prescription given by Home (1986). We extracted data from the 
pipeline science-flatfielded (.fit) 2D files, and used the wavelength arrays of the pipeline extracted (.xld file) 
spectra for calibration. Error arrays were also generated, using the prescription detailed by Home (1986). 
Sensitivity functions were created by dividing the counts by the flux in the pipeline-calibrated spectra. When 
compared, the final optimally extracted spectra appeared very similar to the pipeline data, except that our 
extraction provided: cleaner subtraction of the geocoronal Lya line; the removal of the correct amount 
of background; and the elimination of more hot pixels. Final co-addition required re-sampling individual 
spectra to a common wavelength scale, then adding each spectrum weighted by the inverse of its variance. 



-6- 



3. Ion Column Densities and Element Abundances 

Below, we discuss in detail all of the absorption lines of interest which are relevant to our analysis of 
the metallicity of the interstellar gas in SBS 1543+593. In deriving ion column densities from the analysis of 
the absorption lines, we adopt wavelengths and /-values for the transitions from the extensive compilation 
by Morton (2003). For the Ni II A1317 oscillator strength, we use a new value determined by Jenkins & 
Tripp (2005). When element abundances are compared with their solar counterparts, we use the recent 
reassessment of the solar abundance scale by Lodders (2003). 

3.1. H I 

As can be seen from Figure 1 , we detect two wide damped Lya absorption lines, from the Milky Way 
disk at z — and from SBS 1543+393 at z — 0.009, respectively, as expected from the lower resolution 
spectrum presented in Bowen et al. (2001b). Unfortunately, the wavelength range covered by the G140M 
grating centered at 1222 A is not wide enough to encompass the amount of QSO continuum on either side 
of the two Lya features which is needed to normalize the data. We therefore used the spectra taken with 
the grating centered at 1272 A, which overlap the shorter wavelength data by ~ 4 A, to better define the 
continuum. The two sets of observations used were obtained within a few days of each other (see Table 1), 
so we would expect no variations in the level of the QSO flux between the two spectra. The adopted 
continuum, shown in Figure 1 along with the lcr error to the fit, was derived using the method described in 
Sembach & Savage (1992) which uses a least-squares fit of Legendre polynomials to intensities in regions 
selected either side of a particular absorption line. This results in three normalized spectra: the spectrum 
normalized by the best fit, and two other spectra normalized by the 'upper' and 'lower' 1 a envelopes. 

To measure the H I column density, N(Yl I), we fitted theoretical Voigt profiles to the data. Details of the 
fitting procedures can be found in Bowen et al. (1995). Briefly, theoretical profiles are generated from initial 
guesses of the column density, N, velocity, v, and Doppler parameter, b, of the absorber. After convolution 
with the instrumental Line Spread Function (LSF), the resulting line profile is compared with the observed 
one; a x 2 minimization approach is then used to deduce the best-fitting values of N, b, and v. The LSF 
used was that appropriate for the G140M grating at 1200 A, 2 although the Lya lines in the spectrum of 
HS 1543+5921 are so wide that the shape of the LSF has no influence on the derivation of ,/V(H I). 

We fitted blended profiles simultaneously to both Lya absorption lines from SBS 1543+593 and from 
the Milky Way disk. Wavelengths covering other narrow absorption lines were masked out. There are two 
sources of error in the resulting values of N and v [The shape of the Lya profile is independent of b for such 
a large N(H I)]. The first arises from errors in the continuum fitting. These can be estimated by refitting 
the spectrum as normalized by the 1 a upper and lower continua. The second source of error in the profile 
fits arises from Poisson noise. To estimate the contribution from this source, we performed a Monte-Carlo 



2 see http://www.stsci.edu/hst/stis/performance/spectral_resolution/ 
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1200 1220 1240 1260 1280 

o 

Observed Wavelength (A) 

Fig. 1. — Merged STIS spectra of HS 1543+5921 taken Aug-2003, with the G140M grating centered at 1222 A 
(25.25 ksec total exposure) and 1272 A (39.55 ksec; the gray triangle denotes the wavelength at which the data were 
joined. Data redward of this were obtained at the same epoch as the blue data but represent spectra from only 8 of the 
1 1 orbits which were eventually used. Hence the S/N of the data redward of the triangle is less than that of the final 
combined spectrum used for analysis of the absorption lines.) The damped H I Lya lines due to absorption from both 
the Milky Way and SBS 1543+489 are indicated. Theoretical Voigt profile fits to these features are overplotted, along 
with the continuum adopted (dashed line) and its ±la deviations (dotted). For the purposes of this figure, the data 
have been rebinned two pixels into one. The label "geo" indicates residual geocoronal Lya emission. 

approach, as outlined in Bowen et al. (1995). We used the theoretical best-fit profile as a starting point 
and added to this synthetic spectrum the same amount of noise present in the data. The lines were then 
refitted, resulting in a new set of values of N, v, and b. We repeated this process 500 times, then examined 
the distributions of N, b and v. These distributions were well approximated by Gaussians, and we used the 
limits which encompass 68% of the values as estimates of our random error in these parameters. 

The result of this analysis yielded logJV(H I)= 20.42 ± 0.04 ± 0.01 at a velocity of 2882 ± 16 km s" 1 , 
where the first error denotes the uncertainty in column density arising from continuum placement, and the 
second error represents that arising purely from Poisson noise. The error in v from continuum placement is 
negligible. These results are summarized in Table 2. By adding these errors in quadrature, we adopt an H I 
column density of log A^(H I)= 20.42 ± 0.04. The central velocity of the H I absorption can be compared to 
the value of 2868 km s -1 derived from 21 cm emission measurements (Bowen et al. 2001a). The difference 
of 14 km s" 1 corresponds to approximately one detector pixel for G140M data and is consistent with the 
known errors of 0.5- 1 .0 pixels in the absolute wavelength scale of STIS MAMA observations (Kim Quijano 
2004). However, the velocity of the metal absorption lines discussed below, 2881 km s" 1 , is very close to 
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the Lya absorption velocity; since the metal lines were recorded at different grating tilts and, in some cases, 
at different epochs, the consistency in their velocities suggests that a value close to 2881 km s -1 does indeed 
represent the absorption system velocity, offset some 13 km s" 1 from the systemic velocity of the galaxy, as 
measured from 21 cm observations. 

3.2. S II 

The S II AA1259.52, 1253.81, 1250.58 triplet is clearly detected at the absorption redshift of SBS 1543+593 
(see Figure 2). We define the velocity of the metal absorption line system to be 2881 km s" 1 based on the 
strength, S/N, and symmetry of the Si II A1260 line (see §3.5). With v fixed to this value, Voigt profile fits 
to the S II triplet give N(S II) = 1 .55 ± 0. 14 ± 0.08 x 10 15 cm" 2 and b = 38.5 ± 2.9 ± 2.7 km s" 1 . The latter 
measurement suggests that the S II lines are resolved. As can be seen from the breakdown in the errors, 
the continuum around the S II lines shows variations on scales of a few A, and the resulting continuum 
errors are significant. Again, by adding the errors in quadrature, we arrive at the adopted column density of 
logN(S II) = 15. 19 ±0.04. 

We can compare this value of N(S II) with that derived using the Apparent Optical Depth (AOD) 
method (Savage & Sembach 1991; Jenkins 1996), which measures column density as a function of velocity, 
N a (v). The AOD method provides an important check on the results from the profile fitting analysis: if 
the lines of the S II triplet contain any unresolved, saturated components, the N a (v) profiles of the stronger 
transitions (higher /-values) would appear depressed compared to the weaker ones. A comparison of N a (v) 
for each of the three members of the S II triplet shows no obvious evidence for saturation, to within the 
errors for each line. The column density for each line, integrated over a conservative velocity range of 
-200 < v < 200 km s" 1 , is given in the last column of Table 2. If we again add the separate errors in 
quadrature, and weight each line by its error, we calculate a weighted average of logAf a (S II) = 15. 19 ±0.06, 
which is identical to the value derived from profile fitting. 

Dividing by the neutral hydrogen column density, we find logiV(S II)/iV(H I) = -5.22 ±0.06. Adopting 
the solar abundance log (S /H) = -4.81 ± 0.04 and taking (S/H) = N{S II)/N(H I), we reach the conclusion 
that 

[S/H] = -0.41 ±0.06 

where we have used the usual logarithmic notation 3 . Hence the sulfur abundance in SBS 1543+593 is 
0.39 ± 0.05 times the solar value. We discuss the validity of the assumption that (S/H) = N(S II)/N(H I) in 
detail in §4. 



3 [X/Y]=log(X/Y)-log(X/Y) 
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Fig. 2. — The top three panels show the S II lines arising in SBS 1543+593. Theoretical Voigt profiles which best fit 
the data are overplotted, generated using values of logA^(S II)= 15.19, b = 38.5 km s" 1 , and v = 2881 km s" 1 . (The 
S II A1259 line is blended with Si II A1260, which is fit with a line profile using the values given in Table 2.) The 
bottom panel shows the expected position of Ni II absorption from the galaxy. Although a weak feature may appear 
to be present, it is not statistically significant. 



Table 2. Measurements of absorption lines arising from SBS 1543+593 



Profile fits' 





A 




v© 


b 


\ogN 


\ogN a a > b 


Ion 


(A) 


/-value 


(km s _1 ) 


(km s _1 ) 


(log [cm" 2 ]) 


(log [cm" 2 ]) 


1 1 1 


1 91 S 670 




7CR9-I-1 f. 
ZooZXID 




90 A9-I-0 OzL -1- 01 




S II 


1259.518 


0.0166 








15.18±0.05±0.04 




1253.805 


0.0109 


2881 ±2 


38.5±3.0±2.6 


15.19±0.04±0.02 


15.25±0.10±0.05 




1250.578 


0.0054 








15.18±0.26±0.11 


Ni II 


1317.217 


0.0571 








< 13.83 c 


01 


1302.168 


0.0486 






> 16.2 d 


> 15.3 e 


Si II 


1260.422 
1304.370 


1.1800 
0.0863 




< 41 


> 15.1 


> 15.0 e 


Si III 


1206.500 


1.6300 






> 14.5 d 


> 13.9 e 


CI 


1277.245 


0.0853 








< 13.56 c 


CIT 


1335.708 


0.1283 






14.28±0.15 




Si II* 


1264.738 


1.0500 








< 12.55 c 



a Where more than one error is given, the first error indicates the uncertainty arising from continuum errors, 
while the second indicates that from Poisson noise. 

b Column densities derived from the AOD method, with N a (v) integrated over -200 < v < 200 km s" 1 , where 
v = 2881kms _1 . 

c 2<7 AOD column density limits, derived from continuum and noise errors calculated over a velocity interval 
-80 < v < 80 km s" 1 , where v = 2881 km s" 1 . 

d This lower limit comes from a fit to the line assuming b and v equal that of the S II line given in row 3. 

e Pixels with optical depths r > 3 are set to r = 3 when calculating this limit. 
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3.3. Ni II 

The wavelength at which Ni II A1317.217 absorption is expected from SBS 1543+593 (1329.85 A) 
lies in a relatively clean part of the spectrum. A weak feature may be present at the correct wavelength 
(see Figure 2), but the line is not statistically significant. We measured the AOD column density over a 
velocity interval of -80 < v < 80 km s -1 either side of the velocity at which the Ni II line is expected. 
This velocity range corresponds approximately to the total width of the S II A1250.58 line (see Figure 2). 
Since this S II transition is usually much stronger than the Ni II A1317 line, it is highly unlikely that a Ni II 
line would be wider than the observed S II line. The measured column density over this velocity interval 
is 3.23 x 10 13 cm" 2 . In comparison, the continuum error is 2.13 x 10 13 cm -2 , while the error from the 
noise is 2.64 x 10 13 cm" 2 . Adding these errors in quadrature yields a total 2a error of 6.75 x 10 13 cm -2 , or 
log./V(Ni II) < 13.83. Since the measured column density is roughly half this value, we confirm that Ni II is 
undetected at the 2a level. 

Dividing by N(H I) as before, we find logA^(Ni II)/Af(H I) < -6.59. Subtracting the solar abundance 
log(Ni/H) = -5.78 ± 0.03, and taking (Ni/H) = N(Ni ll)/N(U I), we obtain 

[Ni/H] < -0.81 

We conclude that, in the interstellar gas of SBS 1543+593, nickel is less abundant than sulfur by at 
least a factor of 2.5. 



3.4. O I & N I 

The analysis of the O I A 1302. 2 line required some additional data processing. The spectrum showed 
evidence of defective pixels or another source of non-random noise in the core of the absorption line. Since 
we shifted the target along the long axis of the slit between exposures, we were able to investigate this 
problem by comparing individual spectra lying at different positions on the detector. Of the eight sub- 
exposures of HS 1543+5921 acquired, spectra were recorded at one of three possible locations. We found 
clear evidence for two hot pixels at wavelengths close to the core of the O I line in two sub-exposures 
which were recorded at the same position on the detector (Archive root names o8mr05030 and o8mr06030). 
A third sub-exposure (o8mr05050) also appeared to be unusually noisy (at all wavelengths) compared to 
others. Excluding these three sub-exposures from the total improved the S/N ratio in the region of the 
OI A 1302 line. 

As expected, the O I A1302 line is strongly saturated (see Figure 3). We used both the AOD method and 
Voigt profile fitting to set lower limits to the corresponding column density of O I. In the first case, we set all 
pixels within the line with optical depth r > 3 to r = 3, the limit beyond which we can no longer measure the 
optical depth reliably. Integrating over the same velocity interval as the S II lines, -200 < v < 200 km s -1 , 
we obtain \ogN a {0 I)> 15.3 . 

The second method involves finding the lowest column density which will produce an acceptable fit 
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-500 500 

Velocity relative to 2881 km s" 1 

Fig. 3. — Absorption lines arising from SBS 1543+593 for which, at best, only limits to ion column densities can 
be reliably derived. Top: The O I A 1302 line profile plotted over the data represents a best fit for a component with 
b(0 I)= b(S II) = 38.5 km s -1 . Bottom: With the damped Lya profile from the Milky Way removed, the N I triplet 
has too low a S/N for a reliable column density to be measured. We overplot two sets of profiles, each with the 
same b value as S II, corresponding to the range of column densities of N I implied by the H II region abundance 
determinations of Schulte-Ladbeck et al. (2004): logJV(N I) = 14.9 and 15.4. 

to the observed line profile. Keeping v fixed to that of S II, and allowing b and N to vary, we find that 
b = 49.7 km s -1 and logAffO I)= 15.59. However, there exists one line of reasoning which suggests N(0 I) 
is larger than this. We expect O I to follow H I closely, since charge exchange locks the ionization of 
oxygen to that of hydrogen (Field & Steigman 1971). On the other hand, while we expect most of the 
S II to also arise in the H I material, it is possible that some of the S II might arise in H II gas. The 
velocities of these ionized components may or may not be different from those of O I components, but 
however they are arranged in velocity, we would expect the O I line to be as narrow, or narrower, than the 
S II line, i.e. b(0 I) < b(S II). (For this to be true, [S/O] must be ~ in the absorbing regions, which is 
a reasonable assumption considering that both species are a-capture elements.) This condition is not met 
when N and b are allowed to vary for a fit to the O I line: as stated above, b(0 I) = 49.7 km s" 1 , whereas 
we measured b(S II) = 38.5 km s" 1 in §3.2. We therefore consider that a better limit is obtained by fixing 
b(0 I) < b(S II) = 38.5 ±4.0 km s" 1 . With this constraint, we find that logA^(0 I) > 16.2. This fit is shown 
in Figure 3. If b is uncertain by ±4 km s _1 , so that b(0 I)= [42.5,34.5] km s _1 , then logiV(0 1) = [15.9,16.6]. 
This value of N(0 I) remains a lower limit since b(0 I) could in principle be narrower than b(S II). 
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Adopting logAr(0 I) > 16.2 as the lower limit, we derive logN(0 l)/N(U I) > -4.2. Since (0/H)= 
N{0 I)/N(H I), and log(0/H) Q = -3.31 ±0.05, we find that the lower limit to the oxygen abundance in 
SBS 1543+593 is 

[O/H] >-0.9 

Although we believe the criterion b(0 I) < b(S II) is a reasonable assumption, we should note that 
such a posit is not completely robust. Since the O I A1302 transition is so strong, it is always possible that 
weak O I high-velocity components could be present which have no counterparts in the S II lines, simply 
because the S II /-values are not high enough to produce any absorption. These O I components could blend 
with the main profile in such as way that the measured b value becomes greater than the b value of S II. 
Nevertheless, as can be seen in Figure 2, the O I profile is remarkably close to the theoretical line profile for 
a single component. For high velocity O I components to widen the observed line profile, they would have 
to blend symmetrically on both sides of the main profile, and decrease in strength as their velocity increased 
away from the main absorption complex. 

Our data provide even fewer constraints on the N I abundance in SBS 1543+593. The bottom panel of 
Figure 3 shows the N I triplet after removal of the Milky Way damped Lya profile. The S/N is extremely 
low, and although absorption is present, the lines of the triplet are ill-defined. Schulte-Ladbeck et al. (2004) 
measured log(0/H)= -3.8 ±0.2 from the emission lines of the brightest H II region in the south-west outer 
spiral arm of the galaxy. They also determined that log(N/0) lay between -1.7 and -1.2. It follows that 
log (N/H) lies between -5.5 and -5.0 for gas in the H II region. We can use the value of N(H I) measured 
above, and predict that, if the interstellar gas along the sightline had the same metallicity as the H II region, 
we should measure logN(N I) between 14.9 and 15.4 [again assuming that (N/H) = N(N l)/N(R I)]. Model 
N I lines, given these column densities and the same b value found for the S II triplet, are shown plotted 
over the data in Figure 3. We can say that the data are roughly consistent with these column densities, but 
we cannot distinguish between the upper and lower values of log N(N I). 



3.5. Si II & Si III 

A lower limit to the Si II column density can be derived in a similar way to that described for O I in 
§3.4. Both Si II A1260.4 and 1304.4 lines are detected in the ISM of SBS 1543+593: the 1304.4 A line has 
the smaller /-value and is the least saturated of the two, so should give a more reliable column density; the 
1260.4 A line is more saturated, but its inclusion in the process of fitting a theoretical line profile ensures 
that the b value does not become too wide. With b and N allowed to vary, we find logAf(Si II) > 15.1 and 
KSiH)< 41 kms" 1 . 

We also note the detection of Si III in absorption from SBS 1543+593 which appears in the red wing of 
the damped Lya absorption profile from the Milky Way. The line is recorded at low S/N, but is of some use 
in providing information on the ionization state of the gas, discussed in §4.1 below. Using the AOD method 
to set a conservative lower limit, we arrive at logA^Si III) > 13.9, after again setting all points with optical 
depth r>3tober = 3.A Voigt profile fit to the line keeping v fixed at 2881 km s" 1 and letting b and N 
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Fig. 4. — The C II* A1335.7 absorption line (black,solid) overplotted on the S II A1259.5 line (red,solid), both of 
which arise in the ISM of SBS 1543+593. Although the C II* line falls right at the edge of our spectrum, the portion 
which remains suggests that the line is of similar strength and shape as the S II line. (Other interstellar lines from the 
ISM of SBS 1543+593 are labelled.) A theoretical C II* line profile with the same b and v derived for the S II triplet, 
generated with logA^(C II*) = 14.28 (see Table 2) is shown as a black dotted line. 

vary produces the same Si III column density. We can also set Z?(Si III) < b(S II) as we did for O I; we then 
find logAf(Si III) > 14.5. Although the data are of too low a quality to provide a rigorous estimate of the 
Si III column density, the similarity of these results suggests that our derived lower limit to N(Si III) may be 
robust. 



3.6. C I 

As we discuss in §5, even though no C I A1277 is detected from SBS 1543+593, a limit on the C I 
column density is of interest for deriving limits on the electron density within the H I absorbing gas. The 2a 
column density limit, again measured over the velocity interval of -80 < v < 80 km s" 1 , and combining 
noise and continuum errors in quadrature, is 3.64 x 10 cm , or log/Y(C I) < 13.56; the measured column 
density in the same interval is only 1.58 x 10 13 cm" 2 . 



3.7. C II* 

The C II* A1335 line arising from SBS 1543+593 lies at the very edge of the spectrum taken with 
the G140M grating centered at 1321 A. Unfortunately, only half of the line is recorded (see Figure 4). 
Nevertheless, C II* is extremely valuable in providing insights on the conditions inside the gas clouds which 
give rise to the S II and C II absorption lines. 
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Figure 4 shows that the portion of the C II* absorption we observed has a strength and shape that 
are virtually identical to the corresponding part of the S II A1259.5 line. If the velocity structures of the 
two species are the same, then so are the optical depths of each line, and we can state that N(C II*) = 
(/A)i259.5Af(S II)/(/A) 1335.7 = 10 14 28 , even if the two lines are somewhat saturated 4 . A Voigt profile with 
this column density, and v and b fixed to that of the S II lines, fits the visible line profile quite accurately, as 
can be seen in Fig 4. We can estimate an approximate error on N(C II*) by seeing how the fit changes as N 
changes, again keeping v and b fixed to the S II values. Although this is somewhat more uncertain than a fit 
to an entire line profile, we estimate that N(C II*) is accurate to within ±0.15 dex. 



4. Physical conditions in the ISM of SBS 1543+593 

4.1. Ionization corrections to the neutral gas: CLOUDY modeling 

In the previous sections, we assumed that ionization corrections can be neglected in the derivation of 
abundances. To test whether this is a reasonable assumption, we have assembled photoionization models as 
described in Tripp et al. (2003) using CLOUDY (v96; see Ferland et al. 1998 for a general description of 
the code) with the H I fixed at log N(K I) = 20.42 and [M/H] = -0.41. Since the sight line obviously passes 
through the disk of SBS 1543+593 and starlight is readily apparent nearby (see Figure 2 in Bowen et al. 
2001b or Figure 1 in Schulte-Ladbeck et al. 2004), we have assumed that the gas is photoionized mainly 
by starlight in the disk of a galaxy. To approximate the ionizing flux, we have assumed that the shape of 
the radiation field is similar to that in the diffuse ISM in the disk of the Milky Way, and we have employed 
the Galactic disk radiation field recently presented by Fox et al. (2005, see the "0 kpc" curve shown in their 
Figure 8). Given the large number of H II regions distributed across the disk of SBS 1543+593, this seems 
like a suitable model. It is difficult to estimate the intensity for normalization of the ionizing flux field, 
but since the predicted columns are primarily dependent on the ionization parameter U (= ionizing photon 
density/total H number density = « 7 /«h), the main effect of changing the normalizing flux intensity is to 
change the corresponding «h implied by the model; the relative column densities remain the same when n 7 
and «h are scaled up or down by the same amount. 

Figure 5 compares the photoionization model column densities to observational constraints including 
the measured N(S II) and the 2a upper limit on N(Ni II). We also show the Si III column obtained by fitting 
the Si III profile with the b-value allowed to freely vary (lower value) and with b forced to have the same 
value as the S II lines (upper value in Figure 5). Even the higher value for N(Si III) could significantly 
underestimate the true column due to saturation and an uncertain &-value. 

Figure 5 shows two interesting results. First, even though there can be substantial amounts of ionized 
gas in a gas slab with the SBS 1543+593 H I column that is photoionized by starlight (as revealed by the 
Si III curve), the ionization corrections for S II and Ni II are likely to remain rather small or negligible. The 



4 For the C II* absorption, we assume that the 1335.663 and 1335. 708A transitions are equivalent to a single transition with 
log/A = 2.234. 
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Fig. 5. — Column densities of S II (filled circles), Ni II (open diamonds), and Si III (open squares) predicted by 
a CLOUDY photoionization model for a gas slab with log N(H I) = 20.42, overall metallicity [M/H] = -0.41, and 
solar relative abundances of S, Ni, and Si. The ionizing radiation field was assumed to be similar to the radiation 
field estimated by Fox et al. (2005) for the diffuse ISM in the disk of the Milky Way (see the "0 kpc" curve shown 
in Figure 8 of Fox et al.). The model column densities are shown with small symbols and are plotted as a function 
of the ionization parameter U. The observed column densities are shown with larger symbols at log U « -3.45; the 
hatched region indicates the range of log N(S II) consistent with the observed S II column within ±2a. Two values for 
N(Si III) are shown corresponding to different assumptions about the Si III b-value; conservatively, these should both 
be treated as lower limits (see text). 

reason for this is that while the ion fractions / (averaged over the thickness of the gas slab) of H I, S II, 
and Ni II steadily decrease as U increases, over the entire range of U shown in Figure 5, /(H I) rs f(S II) 
« /(Ni II), so that ionization corrections are not necessary (see eqn. 2 in Tripp et al. 2003). Indeed, the 
strong Si III absorption detected in the SBS 1543+593 DLA is entirely consistent with a simple photoionized 
gas slab, i.e., there is no need based on the current data to invoke a more complex multiphase model. The 
Si III arises in the ionized surface layer of the gas slab, but most of the H I, S II, and Ni II is located in the 
shielded interior of the slab. Second, the underabundance of Ni II is probably not an ionization effect. Over 
a wide range of U, we find little change in N(S II)/iV(Ni II), and at all values of U that we have considered, 
the model predicts a significantly higher Ni II column than the 2a upper limit. 

We must note, however, that this model is not unique. It is straightforward to construct a UV radiation 
field which is much softer than that described above; for example, one in which the flux field is dominated 
by photons from an ensemble of cooler, B-type stars. In this model, hydrogen is ionized over a longer path 
length than for the harder flux field, and some S II arises in H II gas. Whether this UV field is applicable for 
SBS 1543+593 is unclear, however; since over 30 candidate H II regions are identified by Schulte-Ladbeck 
et al. (2004) in the the galaxy, it seems likely that some contribution to the background must come from hot 
O-type stars. Nevertheless, it is a concern that the sightline to the background QSO passes near the center 
of the galaxy, somewhat distant from the most intense areas of star formation. 

A final caveat centers on the fact that it may simply be inappropriate to use a model which assumes that 
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a gas slab is illuminated by a background radiation field, which is the calculation that CLOUDY performs. In 
particular, the state of the ISM illuminated from stars and H II regions embedded in a disk may be different 
from that calculated by CLOUDY. 

Given these uncertainties, in the following section we develop an alternative method for estimating how 
much S II may actually arise in ionized gas. 



4.2. An Alternate Method for Determining the Ionization Correction 

Our ultimate objective is to determine the fraction 

FHn(SII) -^(SII) HI+ iV(SII)Hn (1) 

so that we can apply a multiplicative correction 1 -Fh n(S II) to our measured value of N(S II) and thereby 
determine the correct abundance of sulfur [S/H] in SBS 1543+593. One way to do this is to compare two 
different methods for deriving the expected emission measure (EM) of the ionized gas f n(e) 2 dl in the H II 
regions that are created by radiation from young stars during active star formation. The first procedure, 
discussed in §4.2.1, starts with an estimate for the star formation rate (SFR) per unit area, Xsfr, based 
on our measured N(R I) combined with a widely applicable empirical relationship between Ssfr and gas 
surface density. On the assumption that little if any of the ionizing radiation from the stars is lost, we can 
can convert £sfr into an expectation for the EM (§4.2.2). 

The other, independent method for deriving the EM is through our observation of C IP(§3.7). How- 
ever, the formulation we develop for the EM, derived in §4.2.3, involves a number of previously unknown 
quantities, such as the amount of multiply-ionized C in the H II regions, the apportionment of C II* between 
H I and H II regions, and, finally, the quantity we wish to determine, 1— Fh ii(S II) . By studying the 
production of C II* in H I regions of our own Galaxy, we conclude in §4.2.4 that only a small fraction of 
the C II* in SBS 1543+593 should come from H I regions, after we compensate for differences in the star 
formation rates and dust grain contents in the two systems. Making use of this information, we equate the 
two different determinations of the EM in §4.2.5 and find that we can set a useful upper limit for the amount 
of S II that arises from H II regions. 



4.2.1. Star Formation Rate Based on N( H I) 

Unfortunately, it is difficult for us to determine the SFR near the line of sight by traditional means 
(e.g., stellar brightness and spectral energy distribution, the strength of Ha emission, or the far infrared 
brightness) because we are blinded by the flux from the quasar. Instead, we make use of a recent refinement 
(Kennicutt 1998a) of the Schmidt law (Schmidt 1959) that expresses an empirical relationship between the 
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SFR per unit area Ssfr and the surface density of (atomic plus molecular) hydrogen S gas 

/ £ \ 1 - 4±015 

E SFR = 2.5 ± 0.7 x 10" 4 gas M yr" 1 kpc" 2 (2) 

VlM pc J 

that seems to apply to a broad range of galaxy types and locations within galaxies (including nuclear star- 
bursts associated with molecular gas disks). We see no reason to believe that there could be a breakdown 
of this relationship caused by either prominent bars or activity related to galaxy mergers or tidal effects. 
Drawing from experience in our own Galaxy (Savage et al. 1977) and the Magellanic Clouds (Tumlinson 
et al. 2002), we expect that for our observed value of N(H I) the contribution from molecular hydrogen is 
unlikely to exceed a few percent of the total hydrogen. Thus, for N(H I) = 1 o 20 - 42±0 - 04 cm -2 , we obtain 

Ssfr = 7^ x 10" 4 M yr" 1 kpc" 2 , (3) 

where the stated limits include the effects from the two errors shown in Eq. 2 (which are assumed to be 
uncorrelated), a natural rms dispersion of ±0.3 dex of real outcomes on either side of the Schmidt law, and 
our uncertainty in N(Yl I), all combined in quadrature. 



4.2.2. Anticipated Emission Measure 

Assuming that there is no escape of Lyman Limit photons and that they are not absorbed by dust, 
we can calculate the value of the EM from Esfr- For a Salpeter (1955) IMF and solar abundances, the 
hydrogen-ionizing photon production rate <2(H°) in units of s" 1 is given by 

<2(H°) = SFR/1.08 x 10 -53 M yr" 1 s (4) 

(Kennicutt 1998b). Certain factors could change the constant in Eq. 4. For instance, an IMF defined by 
Scalo (1986) would increase the constant by about a factor of three (Kennicutt 1998b). We acknowledge 
that the stars in SBS 1543+593 will have compositions below solar, but the change in ionizing fluxes should 
be very small (Lanz & Hubeny 2003). 

If all of the ionizing photons are consumed by the ionization of H, <2(H°) per unit area must balance the 
product of EM and the hydrogen recombination coefficient c^ 2) to all levels higher than the ground electronic 
state. For a temperature in the vicinity of 10 4 K (or T4 = r/10 4 K « 1) a power-law fit to the values of a (2) 
and temperature given by Spitzer (1978) yields 

a (2) = 2.55 x 10 -13 r 4 -° 833 cm 3 s" 1 . (5) 

Thus, we find that 

EM = 1 .23 x 10 4 r 4 ° n 8 ( 33 2 S SF R cm" 6 pc (6) 

where, as in Eq. 3, Ssfr is expressed in the units M yr _1 kpc~ 2 , and for an inhomogeneous H II region 
r 4 „ (e) 2 is the temperature (in units of 10 4 K) weighted by the local squares of the electron densities along the 

line of sight, i.e., T 4 n(e) i = J T^n{e) 2 dl I J n(e) 2 dl. 
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4.2.3. C II* Contributions from H II Regions 

Collisions with electrons excite singly ionized carbon atoms to the upper fine-structure level of the 
ground electronic state, creating C II* . The balance between this excitation and radiative decay leads to an 
equilibrium equation from which one can derive the electron density 

g[ A 2d T Q - 5 exp(91K/D { n(CU*) \ 05 .... /n(Cr)\ _ 3 

= 8.63 x iF^ ( KCTT7 J = 183 r 4 ■ exp(0.0091/74) j cm , (7) 

where the lower level's statistical weight g\ = 2, the upper level's spontaneous decay rate A2.1 = 2.29 x 
10" 6 s _1 (Nussbaumer & Storey 1981), and the collision strength tt h2 = 2.90 at T = 10 4 K (Hayes & Nuss- 
baumer 1984). This equation arises from a simplified form of the equilibrium equation in the limit where 
collisional de-excitations can be neglected, which happens when N(C II*)/iV(C II) <C 1 (which is true here) 
- see, e.g., Jenkins et al. 2000. Henceforth, we will drop the exp(0.0091/r4) term since it is so close to 1.0 
when T4 ~ 1, which seems justified for H II regions in systems with metal contents somewhat below that of 
our Galaxy, such as those in M101 (Kennicutt et al. 2003). In most circumstances, optical pumping of the 
excited level is unimportant (Spitzer & Jenkins 1975; Sarazin et al. 1979). 

An alternate expression for n(e) within the H II region is given by the equation 

nW=Qn(cn)j(cn) (8) 

where 

rc(cn)+n(cm)+n(crv) 
y(C ii) = -77777 (9) 

n(C II) 

In an H II region where the ionization by starlight dominates over that from an intergalactic field, the 
expected values of y(C II) are strongly dependent on the average temperature of the exciting stars : models 
of ion abundances in H II regions evaluated by Stasinska (1990) indicate that the ratio of 
/ n(C total )n(e)dV/ J n(C U)n(e) dV varies from 1.26 for T* = 32,500 K, to 8.5 for T* = 40,000 K, and could 
reach as high as 22 for T* = 50, 000 K. 

It is reasonable to assume that [S/C] « in SBS 1543+593, which permits us to arrive at an estimate 
for the H to C abundance ratio, 



cj Vc /0 



l-F Hn (S II) 



N(S II) 



(10) 



which we consider to be a universal abundance ratio everywhere. We can now formulate a useful expression 
for 72(e) 2 by taking the product of the expressions in Eqs. 7 and 8, 

EM = 5.93 x 10 -18 J n(C II*)r 4 - 5 v(C 11) dl cnf 6 pc (11) 

To be rigorous, we must allow for the possibility that T4 and y(C II) could vary inside the H II region(s), and 
therefore we define a product 

Yn(c n*)(C II) r 4 ° n 5 (C ir) = j T 4 5 y(C II) n(C II*) dl J n(C II*) dl (12) 
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so that we can work with a more streamlined representation 

fN(U I) 



EM = 5.93X10-I - vw(sn) 



Y n (C n*)(C II) F 4 ° n 5 (C n *)Fn n(C II*) 
1-WS ID 



JV(C II*) croupe (13) 



after we make the substitution for (H/C) given in Eq. 10. [Fh n(C II*) is defined in a manner identical to that 
for S II in Eq. 1.] 

Were it not for the uncertainties in the quantities y„ (C n*)(C II), T®^ c u *y and Fh n(C II*), we could 
determine Fh h(S II) by equating the EM given in Eq. 13 with that given earlier in Eq. 6. To provide useful 
constraints on these unknown quantities, we must evaluate the probable contribution of C II* from H I 
regions. 



4.2.4. An Evaluation ofF u n(C II*) 



The discussion presented in the previous subsection focused on the properties of fully ionized gas 
regions, as revealed by C II*, but with the retention of an unknown fractional quantity Fh n(C II*) that 
expressed how much of the observed C II* actually came from such regions. We now examine the plausi- 
bility that Fh ii (C II*) could be appreciably lower than 1.0, and that a reasonable fraction 1 -Fh h(C II*) 
of the C II* arises from H I regions, as opposed to H II regions. One way to do this in principle is to 
measure N(Si H*)/N(C II*) (Howk et al. 2005), since the excitation energies of the two are very differ- 
ent. Unfortunately, our 2a column density limit derived at the expected position of the Si II* A 1264 line is 
logAf(Si II*) < 12.55. This limit is not sufficiently sensitive to discriminate between contributions from H II 
regions with T4 ~ 1 and H I regions at much lower temperatures. 

We turn now to a different, but less direct way to measure Fh h(C II*). For the H I regions that intersect 
our line of sight, the average cooling rate per H atom arising from the emission of radiation at 157.7 fim 
caused by the radiative decay of C II* corresponds to 



k = (A 2 ,ihc/X) 



= 2. 



x 10 
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l-F Hn (C II*) 
1-Fhh(C II*) 



N(C H*)/iV(H I) 
N(C n*)/JV(H I) erg s" 1 H atom -1 . 



(14) 



For our measurements of N(C II*) and N{H I), we obtain log I lc/ 



l-F Hn (C II*) 



) . -25. 



68. If we omit 



the factor 



l-F Hn (C II*) 



, this value is about the same as the rate log/c = -25.70 (+0.19, -0.35 for the 
la dispersion of individual determinations) found by Lehner et al. (2004) for low- velocity clouds and some 
accompanying ionized material in the local region of our Galaxy. 

For densities «(H) > 0.1 cm -3 , within the ISM of our Galaxy the heating rate from energetic photo- 
electrons liberated from dust grains dominates over those arising from cosmic ray and X-ray ionizations 
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(Wolfire et al. 1995; Weingartner & Draine 2001b). To a first order of approximation, we expect the dust 
grain heating to scale in proportion to both (1) the local starlight density, with special emphasis on the far- 
UV photons, and (2) the density of dust grains per H atom. This relationship is valid as long as there is no 
appreciable loss of efficiency caused by the grains acquiring a positive charge because the recombination 
rate with free electrons is too low - see Fig. 16 of Weingartner & Draine (2001b). A reasonably good gauge 
of the local far-UV starlight flux is the SFR. We can propose that the grain density should scale in proportion 
to the overall metal abundance, 10 [M//H] , although we caution that there are indications from IR emission 
that the character of dust in low metallicity systems, and in particular the polycyclic aromatic hydrocarbons 
(PAHs) which are important for heating the gas, is different from that of our Galaxy (Draine 2005; Engel- 
bracht 2005). Thus, there may not be a simple downward scaling of the abundances of all the different kinds 
of dust. Nevertheless, if we overlook this possible complication and use our local region of the Milky Way 
(MW) as a comparison standard, we can express a proportionality 



Z C (SBS 1543+593) = [M/H] Z C (MW) 
S SF r(SBS 1543 + 593) S SF r(MW) ' 



(15) 



where Z C (SBS 1543 + 593) is defined in Eq. 14, S SF r(SBS 1543 + 593) is calculated in Eqs. 6 and 13, 



and 10^/H] = 



1- Fuu(S II) 



10 41 (see §4.2.2). For / C (MW) we adopt a representative value of 1.3 x 

10~ 26 erg s -1 H atom -1 , which Lehner et al. (2004) found for low velocity clouds that have sufficient column 
densities to make the relative contributions of C II* from ionized regions relatively small (see their Fig. 8). 
An acceptable range for Esfr(MW) in our region of the Galaxy is (3.5-5) x 10~ 3 M Q yr" 1 kpc~ 2 (Rana 
1991). Using these quantities, we solve Eq. 15 to find that 



WC II*) = 



i +0.057r-°§ 3 r 4 ° n 5 (C in Y n(C ir)( c id 



(16) 



With some loss of rigor for nonuniform regions where both T4 and y(C II) might vary, we can consolidate 
the two temperature terms with the different weightings into a single representation T4 to obtain 

-1 



F Hn (C II*) = 



l + 0.057r 4 - 0333 F„ (C m(C II) 



(17) 



4.2.5. An Evaluation ofY n( Q n*)(C II) 

As we stated earlier (§4.2.3), the value of the inverse of the fraction of C atoms in the singly ionized 
form, Y n ( C n*)(C II), can be slightly greater than 1 or very large, depending on the hardness of the ionizing 
radiation (i.e., the characteristic temperature of the exciting stars). Any constraints that we can place on the 
value of this parameter are critical for our being able to estimate how much of the C II* arises from H II 
regions (see Eq. 17 above), and this quantity is essential for our final estimate of Fh n(S II). 

If we convert the Ssfr derived from the Schmidt law (Eq. 3) to an EM using Eq. 6 and equate the result 
with the EM given in Eq. 13 (with the expression in Eq. 17 substituted for Fh n(C II*)), we find that 



Y n(C H*)(C II) _ +i.o T 0.333 



l+0.057r-°- 333 F„ (C i n (C II) 



l-F Hn (S II) 



(18) 
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(once again we have consolidated the slightly different kinds of temperature terms T 4 ^ and T^c n*) 
into a single term T4). The uncertainty for the numerical coefficient on the right-hand side of the equation 
considers the uncertainty in the value of Esfr combined with the uncertainties in the column densities of 
H I, S II and C II* that also appear in this result. Clearly, with T4 = 1 we are unable to find an acceptable 
result with the preferred value of the coefficient, since we require that Y n ( C n*)(C II) > 1 .0 and F H n(S II) > 0. 
If this coefficient is at its (la) upper limit of 1.7 and Y n( c n*)(C II) = 1.0, then Fh n(S II) could equal 0.44. 
However, Y n ^ c n*)(C II) is likely to be somewhat greater than 1.0, but by an amount that is uncertain, so 
Fh n(S II) is probably lower than 0.44. 

To summarize, our consideration of the Schmidt law in conjunction with our measurement of N(H I) 
leads to an implied SFR which, when combined with our interpretation on the origin of C II*, places an 
upper limit of about 44% on the relative amount of S II that could arise from H II regions. 

While one might argue that in principle Ssfr could deviate above the +la limit that was used to derive 
our upper limit for Fhh(S II) if the conversion rate of gas into stars were significantly higher than implied 
by the scaling of eq.(2), this possibility seems unlikely for a galaxy with a surface brightness as low as 
that of SBS 1543+593. The argument presented here, along with the one given in §4.1, means that [S/H] 
for the gas we detected in SBS 1543+593 is probably very close to the more naive determination based on 
log[N(S ll)/N(U I)]-log(S/H) Q = -0.41, but it might be as low as -0.66. 



5. Physical Conditions of the Gas 

If we knew the column density of C II in the H II region, we could determine a representative electron 
density weighted by the density of C II inside it, i.e., n(e) n (c n> This would be done through the use of 
Eq. 7 (but with a replacement of local densities n with column densities N), because we know that most of 
the C II* arises from collisions with electrons within fully (or mostly) ionized material. Unfortunately, the 
information about Fh h(S II) cannot be applied to C, because the low value of Fh n(S II) could have arisen 
from the fact that y(S II) is very much greater than 1 over most parts of the region. Indeed, the models of 
Stasinska (1990) indicate that S reverts to multiply ionized stages more strongly than C, and this could in 
principle help to explain the small value of Fh h(S II). 

For the neutral material, the situation is different. Since it is virtually certain that most of the C and S 
is singly ionized, we can make the reasonable assumption that [C/S] = and state that inside the H I region 



wn)-f§) 







l-F Hn (S II) 



N(S II) (19) 



or 

log/V(C II) = 16.39 if F H n (S II) » (20) 

In §3.6, we reported a limit log/V(C I) < 13.56. We can use this limit together with our inferred value for 
log/V(C II) to place limits on n(e) and «(H) in cool H I regions if we use the equilibrium equation for the 
balance of photoionizations, at a rate Tc, against recombinations with free electrons (with a rate constant 
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Fig. 6. — Upper bounds that satisfy Eq. 21 and our observation that logA^(C I)-logA^(C II) < -2.83, i.e., combinations 
of n(e) and n(H) below the curves for four representative temperatures are the only ones allowed by our observations. 
The left-hand panel depicts the condition that only 0.56 times our assumed column density of C II [derived from 
N(S II) and our assumption that [C/S] = 0] is in the H I region, while the right-hand panel shows the conditions when 
virtually all of the C II is in the H I region. The shaded portion in the lower right portion of each panel indicates a 
disallowed region where n(e) is less than the product of n(H) and the sum of the relative abundances of all elements 
that can be photoionized in an H I region. 

a e ) and negatively charged dust grains (with a rate constant a g ), 

[a e n(e) + 1 [M/H] a g n(U)] n(C II) = r c «(C I) . (21) 

Formulae for the rate constants a e and a g are defined by Aldrovandi & Pequignot (1973) and Weingartner 
& Draine (2001a), respectively. As in §4.2.4 above, when we evaluate this equilibrium we assume that 
the factor 10 [M / H] reflects a reasonable estimate for the grain density relative to that in our Galaxy. The 
representative ionization rate for neutral carbon Tc should scale in proportion to Esfr since the atoms are 
ionized by photons with energies greater than 1 1 .26 eV. Thus, once again, we can scale a physical process 
using the Milky Way as an example, 

T C (SBS 1543 + 593) = r c (MW)S SF R(SBS 1543 + 593)/£ SF r(MW) , (22) 

where T C (MW) = 2.24 x lO^V 1 (Jenkins & Shaya 1979). 

Figure 6 shows the upper bounds for the various combinations of n(e) and «(H) that satisfy Eq. 21 
and the condition that the observed logA^C rt-logA'XC II) < -2.83, representing the case Fu n(C II) = 0, 
or the alternative that Fh h(C II) could be as low as our limit 0.44 for Fh h(S II). In effect, these two cases 
recognize that either a modest fraction or virtually all of the observed C II could arise from any H I region 
that could conceivably hold some C I. 
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6. Summary & Discussion 

The primary goal of our STIS program was to obtain the abundance of sulfur and nickel in the in- 
terstellar medium of the galaxy SBS 1543+593, by analyzing ultraviolet absorption lines produced by the 
galaxy in the spectrum of the background QSO HS 1543+5921. With the good quality data obtained, and 
a careful analysis of the absorption lines detected (including an optimal extraction technique implemented 
to improve the standard pipeline extraction), we find that the abundance of sulfur is approximately two- 
fifths solar, [S /H] = -0.41 ± 0.06 (Table 3). This value is in good agreement with the oxygen abundance, 
[O/H] = -0.49 ± 0.2 deduced by Schulte-Ladbeck et al. (2004) and [S /H] ~ 0.3 ± 0.3 measured by Schulte- 
Ladbeck et al. (2005) from an analysis of nebular emission lines in the optical spectrum of the brightest H II 
region in SBS 1543+593. The most straightforward conclusion is that the two methods give a concordant 
picture of the degree of metal enrichment in this galaxy, at least as far as the a-capture elements — which are 
mostly the products of Type II supernovae — are concerned. 

As is often the case, there are caveats to this conclusion. Emission and absorption line measurements 
do not refer to exactly the same location within the galaxy, but the two sight lines are separated by only 
16.6 arcsec, or 3.3 kpc at the redshift of the galaxy, and low luminosity galaxies like SBS 1543+593 
are found to exhibit relatively flat radial abundance gradients (de Blok et al. 1996). In deriving the sulfur 
abundance in the DLA, we have assumed that N(S II)/A^(H I) = (S/H). This assumption would be incorrect 
if sulfur were depleted onto grains — a possibility which we consider far-fetched given the composition of 
local interstellar dust — or if some of the S II absorption arises in ionized gas which does not contribute to 
the damped Lya line. If we have under-estimated the total amount of hydrogen along the sight line, then the 
measured abundance is only an upper limit. 

Testing this last assumption has been a principle aim of this paper. Since O I tracks H I so closely in 
interstellar gas, an accurate measure of [O I/S II] would enable us to measure the fraction of S II arising in 
the observed neutral hydrogen (assuming [O/S] = 0). Unfortunately, we have only a lower limit to N(0 I) 
from the saturated O I A1302 line (§3.4), and can show only that [O I/S II] > -0.49. Hence we can only be 
certain that > 32% of the observed S II line arises in neutral gas. Even the lower limit to N(0 I) is subject 

Table 3. Summary of Abundances in SBS 1543+593 



X 


log(X/H) a 




[X/H] 


s 


-4.81 ±0.04 




-0.41 ±0.06 


o 


-3.31 ±0.05 


> 


-0.9 


Si 


-4.46 ± 0.02 


> 


-0.8 


Ni 


-5.78 ±0.03 


< 


-0.81 



a Reference solar abundances taken from 
the collation by Lodders (2003). 
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to uncertainty: extra components of high velocity O I (relative to the bulk of the absorption) could arrange 
themselves in such as way that we adopt too narrow a Doppler parameter in calculating the lower limit, 
leading us to over-estimate the O I column density. 

We can perhaps look to our own Galaxy to understand how the ionization stages of sulfur track neutral 
hydrogen. In the Milky Way, (S II/H I) ratios through the disk are exactly [S/H] Q (Lehner et al. 2004) at the 
same H I column densities measured towards SBS 1543+593, suggesting that all the S II arises in Milky Way 
H I gas (assuming that most of the ISM does not have super-solar sulfur metallicities). Although conditions 
in the ISM of SBS 1543+593 may be different from that present in the Milky Way, it seems somewhat 
reassuring that at least in one high-Af(H I) environment, most of the S II comes from H I regions. 

Nevertheless, a more quantitative estimate of any required ionization correction would be preferable. 
Hence we have used the CLOUDY photoionization code to estimate how much S II arises in neutral or 
ionized regions. Not surprisingly, we find that the answer depends on the shape of the UV background field 
adopted for the calculation. If a field similar in shape to that of the Milky Way is used, then no ionization 
corrections are required. A much softer background though, or many embedded H II regions, could leave 
some of the S II arising in H II gas. Given the distribution of candidate H II regions in the galaxy identified 
by Schulte-Ladbeck et al. (2004), a Milky Way-like disk model seems plausible, but we acknowledge that 
without a real measure of the UV field near the center of the galaxy, CLOUDY models may only be a 
simplistic approximation of the conditions in the interstellar gas. 

Given all these difficulties, we have developed a different approach to deriving an ionization correction, 
based on comparing two different methods for deriving the emission measure of the ionized gas (§4.2). We 
have used the measured N(H I) to estimate the star formation per unit area in SBS 1543+593, and from 
there — assuming none of the Lyman Limit photons produced by the stars escape (or are absorbed by dust) 
— we derive the emission measure. We have also detected C II* absorption in the ISM of the galaxy, 
which enables us to calculate the same emission measure. Combining these two results yields a relationship 
between the amount of sulfur arising in ionized gas and the amount of carbon which exists as C II. Using 
plausible values for the latter value leads us to believe that no more that 44% of the S II is likely to arise in 
ionized gas. At most, therefore, the sulfur abundance might be as low as -0.66, only 0.25 dex smaller than 
if we had assumed all the S II came from neutral gas. 

Since we do not detect the Ni II A1317 line, we can only deduce an upper limit to the abundance of 
nickel, [Ni/H] < -0.81 or (Ni/H) < 1/7 of solar. This limit implies that nickel is less abundant than sulfur 
by a factor of 3 or more. It seems likely that nickel is underabundant due to depletion by dust, as seen in the 
local ISM, where it is ~ 10-300 times below its solar value (e.g. Savage & Sembach 1996; Jenkins 2004). 
On the other hand, since nickel is an Fe-peak element synthesized mostly by Type la supernovae, it may 
be intrinsically less abundant than sulfur (and oxygen). We might expect there to be little over-abundance 
of alpha elements relative to iron-peak elements, since the relatively high metallicity of sulfur should imply 
that sufficient time has passed for Type la supernovae to build up the iron-peak abundances. Of course, both 
factors may contribute to the underabundance of nickel; in order to assess their relative importance we need 
to measure an undepleted Fe-peak element like zinc. 
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Assuming that our measured sulfur abundance is correct, it is reassuring to confirm that emission- 
and absorption-based abundance measurements do give consistent results, as we would expect. Of course, 
the case reported here refers to only one (late-type dwarf) galaxy but, as emphasized above, the test we 
have carried out is the 'cleanest' so far, free from the complications introduced by: (i) differences in dust 
depletion and nucleosynthetic origin of the elements being compared, (ii) absorption line saturation, and (iii) 
the geometry of composite sight lines through the absorbing medium. Our result, therefore, has implications 
for the interpretation of similar comparisons which have been reported recently, as we now briefly discuss. 

There seems to be no need for an envelope of unprocessed, neutral gas surrounding this dwarf galaxy 
to contribute to the observed absorption, in contrast with the suggestion put forward by Cannon et al. (2005) 
in their interpretation of FUSE absorption line spectroscopy of a number of star-bursting dwarf galaxies. 
Possibly, this difference is related to the evolutionary status of the galaxies involved — SBS 1543+593 is 
more metal-rich than some, though not all, of the galaxies considered by Cannon et al. (2005) and is not 
undergoing a current burst of star formation. One could think of plausible scenarios where the existence of 
a metal-poor halo around galaxies may be related to both of these factors. Of course, measuring abundances 
in dwarf galaxies with FUSE has its own set of unique problems. For one thing, absorption lines are seen 
against the background light of the central starburst which is an extended, rather than point-like, source. 
In addition, the relatively small telescope aperture of FUSE can only record such spectra at modest signal- 
to-noise ratios (S/N); in these conditions the absorption lines studied tend to be strong and saturated. Had 
our metallicity measurements in SBS 1543+593 been based on O I A 1302 and other similarly saturated 
absorption lines, we may well have (erroneously) concluded that there is an order of magnitude difference 
between H I and H II region abundances in this galaxy too (§3.4). 

As for the population of DLAs as a whole, it can be seen from Figure 7 that SBS 1543+593 is among 
the least metal-poor DLAs at z < 1 (and indeed, at any redshift). This may well be the result of the unusually 
low impact parameter of this QSO-galaxy pairing, as the sight line to HS 1543+5921 intersects the galaxy 
at < 0.5 kpc from its center. In other low-z DLAs, radial metallicity gradients have been proposed to 
explain the finding that [Fc/HJdla is lower than [0/H]h ii (Chen et al. 2005). We would have found a similar 
result ourselves, had we based our abundance comparison on nickel (§3.3) which, like iron, can be depleted 
from the gas-phase of the ISM by significant, and variable, amounts. Furthermore, the systematic offsets 
which have recently been identified in the derivation of [O/H] from the ratios of strong nebular emission 
lines (Garnett et al. 2004a,b; Bresolin et al. 2004) seem to be less important when the oxygen abundance 
is sub-solar, as in SBS 1543+5921, than solar or super-solar, as is the case for the galaxies considered by 
Chen et al. (2005). These offsets would work in the sense of increasing the difference between H II and H I 
region abundances. Thus, in order to assess the degree to which abundance gradients are responsible for the 
generally low metallicities of DLAs, we need: (i) better determinations of the gradients in local galaxies and 
(ii) more instances like the one reported here where emission- and absorption-deduced abundances can be 
compared without the uncertainties associated with dust depletion and strong emission line calibrators. 

The cooling rates l c of the DLAs derived from C II* absorption lines have been investigated in detail by 
Wolfe et al. (2003b,a, 2004). The median value of log(/ c ) in the DLA population is ~ -27.00 (where l c is in 
units of ergs s" 1 H atom -1 ). Calculating l c for SBS 1543+593 relies on knowing the fraction of C II* in H II 



-27- 



0.5^ 




—2.5 c i i i i 

12 3 4 

z 

Fig. 7. — The current sample of zinc abundance measurements in DLAs, as a function of redshift z, from the com- 
pilation by Kulkarni et al. (2005). The star close to z = represents the sulfur abundance we have measured in 
SBS 1543+593. Any offset between metallicities measured from zinc and sulfur is likely to be less than 0.2 dex (see 
Figure 9 of Nissen et al. 2004). 

gas, i.e., Fyi n(C II*) in eqn (14). We can use plausible values of Y n (c n*)(C II) to estimate Fh h(C II*) via 
eqn 17: for the range Y n(C ir) (C II) =[1,22] (see §4.2.3), we find log(/ c ) = [-26.9,-25.9]. However, if the 
star formation rate in SBS 1543+593 is close to the value suggested in §4.2.1, ~ 1 -2 x 10" 3 M yr" 1 kpc~ 2 , 
then: Y n ( C n*)(C II) must be between « 1 -2; the fraction of S II which resides in ionized gas must be 
<< 44%; the fraction of C II* arising in ionized gas would be > 90 %; and the cooling rate would be close 
to the median value of log(Z c ) = -27.0 for the DLAs. 

Whether such a comparison is entirely valid is, however, unclear. For the DLAs, C II* absorption is 
thought to come primarily from the cold neutral medium, i.e. with Fh h(C II*) ~ (Wolfe et al. 2004). 
Unfortunately, a more detailed comparison between the methods used to measure l c in this paper and those 
used for the DLA population is beyond the scope of this paper. 

Our observations of SBS 1543+593 provide a strong incentive for finding more examples of close QSO- 
galaxy pairs at low redshifts, where the galaxy producing the DLA in the QSO spectrum can be imaged and 
its properties measured. While on-going large scale surveys, such as the Sloan Digital Sky Survey, have 
dramatically increased the numbers of such pairs, the current lack of a sensitive, high resolution, near-UV 
spectrograph in space will hamper progress in this area until a worthy successor of STIS is flown. 
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